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Abstract 

The role of naturally varying vegetation in influencing the climate variability in the 
Sahel is explored in a coupled atmosphere- land- vegetation model. The Sahel rainfall 
variability is influenced by sea surface temperature (SST) variations in the oceans. 
Land-surface feedback is found to increase this variability both on interannual and 
interdecadal time scales. Interactive vegetation enhances the interdecadal variation 
significantly, but can reduce year to year variability due to a phase lag introduced by 
the relatively slow vegetation adjustment time. Variations in vegetation accompany 
the changes in rainfall, in particular, the multi-decadal drying trend from the 1950 s 
to the 80 s. 

Figure 1A shows the observed rainfall over the West African Sahel region (1) from 1950 to 
1998. It shows a multi-decadal drying trend from the 1950s to the 80s and early 90s, as well as 
strong interannual variability. Proposed mechanisms for this dramatic trend include global sea 
surface temperature variations (2-5) and land-use change, i.e., the desertification process (6-7). 
Since vegetation distribution tends to be largely controlled by climate (8-9), and surface property 
change can impact the climate by modifying the atmospheric energy and water budget (10-12), it 
is reasonable to conjecture that dynamic vegetation-climate interaction might influence the climate 
variability substantially in a climatically sensitive zone like the Sahel. We test this hypothesis in a 
coupled atmosphere-land- vegetation model of intermediate complexity. 

The atmospheric component of the coupled model used here is the Quasi-equilibrium Tropical 
Circulation Model (QTCM) (13, 14), which is coupled to the land-surface model Simple-Land 
(SLand) (14) which represents the first-order effects relevant to climate simulation. The QTCM 
simulates a seasonal climate over the Sahel that is close to observations and that compares favorably 
with current atmospheric general circulation models (GCMs). We model the major effects on cli- 
mate of a varying vegetation through its control of the evapotranspiration process and modification 
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of the surface albedo. Other effects, such as surface roughness and modification of soil properties, 
are not considered here. 

Vegetation growth in the tropics is less influenced by temperature and nutrient limitation, but 
is mainly controlled by water availability on the interannual and decadal time scales because of 
the relatively large rainfall variability. The central equation in the dynamic vegetation model is a 
biomass equation driven by photosynthesis and vegetation loss: 

^ = a/3(w)(l - e~ kL ) - V/t (1) 

where a is a carbon assimilation coefficient, /? is the soil moisture dependence as used in the original 
SLand, L is the plant leaf area index (LAI), and k is the extinction coefficient of photosynthetically 
active sunlight. The vegetation time scale r is taken as one year in the following experiments. This 
is similar to the biomass equations used in models with more explicit vegetation dynamics (15-18). 
V is interpreted as vegetation amount or leaf biomass and it is normalized between 0 and 1. The 
LAI is assumed directly proportional to V as L — L max V , where L max is a maximum LAI taken 
as 8. The current version does not explicitly model plant competition, nor does it consider species- 
specific characteristics such as resources allocation. The seasonality is not explicitly modeled for 
V, so Eq. 1 represents variation on the background of a mean seasonal cycle. 

The original version of the land model is modified to account for the effects of leaf-to-canopy 
scaling (19) such that the canopy conductance for evapotranspiration is: 

9c — 9s mix P{w){l — e~ kL )/k (2) 

where # Smax is a leaf-level maximum conductance. Note that photosynthesis and evapotranspiration 
are closely related in Eq. 1 and Eq. 2 (19). Besides modifying evapotranspiration through Eq. 2, 
vegetation also changes land-surface albedo A as (20) 


A = 0.38 - 0.3(1 - e~ kL ). 


( 3 ) 
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This corresponds to an albedo of 0.38 at V — 0 (desert), and 0.09 at maximum vegetation V = 1 
(dense forest). Thus, vegetation feeds back into the atmosphere by modifying evapotranspiration 
and surface albedo through Eq. 2 and Eq. 3. 

In order to identify the relative importance of oceanic forcing, land-surface and vegetation 
processes, we perform a series of model experiments, starting from a run in which both land and 
vegetation are interactive (the default and realistic case), designated as Lon Von. The monthly 
output from this run is then used to derive a vegetation climatology (with a seasonal cycle but 
which does not change from year to year) which is used as a model boundary condition for the second 
run LonVofF. The output of the run LonVoff is then used to derive a soil moisture climatology which 
is used to drive the third experiment LofFVoff. In all these three runs, the coupled atmosphere-land- 
vegetation model is driven by the observed monthly SST from 1950 to 1998 (21). All three start from 
an initial condition taken from an interactive land/ vegetation run forced by a climatological SST. 
The model simulated annual rainfall over the Sahel from these experiments is shown in Fig. 1B-D. 

Compared to the observations in Fig. 1A, the LofFVoff run with non-interactive soil moisture 
and vegetation (Fig. IB) shows a weak interannual variation, and a much weaker interdecadal 
signal, although a drying trend can be seen from the 1950s to the 80s. The interactive soil moisture 
(Fig. 1C) appears to somewhat increase this interdecadal drying trend all the way into the early 90s. 
The simulated soil moisture shows a high degree of correlation with precipitation. The amplitude 
of interannual variation is also larger in general, in agreement with other studies (22-24). 

By allowing vegetation feedback to the atmosphere, the decadal rainfall variability is enhanced 
significantly (Fig. ID). The wet periods in the 1950s and the dry periods in the 70s and the 80s stand 
out and are much more like in the observations. Compared to the non-interactive vegetation case 
LonVoff, the interannual variability does not show apparent enhancement. In some cases, the year 
to year change such as from the dry 1987 to the relatively wet 1988 is actually reduced, because the 
vegetation is still low from the previous drought due to its relatively slow response. This complicated 
lagged relationship between vegetation and precipitation is also seen in observations (25), but our 
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understanding of vegetation dynamics and the modeling tools available are at present not sufficient 
for a precise assessment. The variation in the rainfall drives a similar trend in vegetation through 
vegetation growth or loss (Fig. ID). The vegetation lags slightly behind the rainfall and its variation 
is also smoother, although these tendencies are not very strong because the one year vegetation 
time scale used in the model runs is comparable to the time resolution of the plot. 

In order to assess the internal climate variability in the model, three 9-member ensemble runs 
have been conducted corresponding to the cases above. The member runs in an ensemble differ 
only in their initial atmospheric and soil moisture conditions. We use the Sahel rainfall difference 
between the wet period 1950-67 and the dry period 1968-87 as an indicator of the amplitude of 
the interdecadal variation. Figure 2 shows the successive amplification of the decadal trend with 
the inclusion of interactive soil moisture and especially, vegetation. However, even the Lon Von 
case still tends to underestimate the observed decadal trend. Furthermore, the scatter among the 
ensemble members also increases when additional feedbacks are included. Interactive vegetation 
increases the variance even though vegetation initial conditions are identical in these runs. 

The interactive vegetation modifies the precipitation through a chain of positive feedback loops. 
For instance, decreased rainfall leads to less water availability and reduces vegetation, which in 
turn leads to higher surface albedo and reduced evapotranspiration. This weakens the large-scale 
atmospheric circulation by reducing the energy and water flux into the atmosphere column, thus 
further decreasing the local rainfall (6, 10, 12, 26). 

The dynamic nature of the vegetation-climate interaction can be understood more precisely in 

a linear system by simplifying these feedback processes: 

dV f aP'-V f , A \ 

~ir = (4) 

at r 

P' = fiV' + F 0 e iwt (5) 

Here V' and P' are perturbations to vegetation and precipitation respectively, and they approximate 
the interannual/interdecadal anomalies shown above. The coefficients a and n represent the 
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strength of the local interaction between vegetation and the atmosphere. Equation 4 is a linearized 
version of the biomass equation (Eq. 1) because soil moisture is forced by precipitation, and Eq. 5 
approximates the vegetation feedback to rainfall through changes in surface albedo in Eq. 3 and 
evapotranspiration in Eq. 2. The rainfall is sinusoidally forced at frequency u and amplitude F 0 , 
representing the SST-induced change in the large-scale atmospheric circulation. 

The dependences on u of the amplitude of P' and phase lag <j) are shown in Fig. 3 (27). At low 
frequency forcing (lot -+ 0), vegetation has enough time to establish a near equilibrium with the 
precipitation, consequently the precipitation is enhanced by a factor of 1/(1 - ay) with little phase 
lag. This explains the amplification of the interdecada! variation of rainfall shown in Fig. ID. At 
high frequency forcing (wr — ► oo), vegetation has little time to respond to the forcing due to its 
relatively long adjustment time, therefore, the precipitation variation is mostly a direct response to 
the forcing. At intermediate frequency forcing, the amplitude is enhanced slightly, but the phase 
lag is at a maximum. This phase lag has significant consequences for the interannual variability. In 
1988, for instance, the memory in vegetation of the previous drought years has delayed and reduced 
the otherwise strong wetting tendency resulting from SST-induced atmospheric circulation change 
(cf. Fig. 1C and D). 

The model simulated Sahel rainfall in Lon Von shows a correlation with the observation of 0.67, 
a significant improvement from a 0.44 correlation in LonVoff. However, the year to year comparison 
with the observations is not as satisfactory. When we decompose the Sahel rainfall time series into 
low frequency (longer than 10 years) and high frequency (shorter than 10 years) components, the 
correlation with the observation is only 0.1 at high frequency and 0.94 at low frequency for the 
Lon Von run. The discrepancy in the interannual simulation, both in our model and in the GCMs 
(4, 28), is not necessarily due to model deficiencies, since considerable variability arises from the 
chaotic nature of the atmosphere-land-vegetation system. 

While we focus on natural climate variations involving vegetation changes, this does not exclude 
any role anthropogenic land-use change might play. Indeed it is possible that the anthropogenic 
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factor might account for the remaining difference between the interactive vegetation run and 
observations in Fig. 2 on decadal time scales. But because the vegetation feedback acts to amplify 
Sahel rainfall variability that originates from SST variations, significant effects can occur with 
relatively small vegetation changes. In the experiment with interactive vegetation (Fig. 3D), the 
change of surface albedo (not shown) on the interdecadal time scale is about 0.03. This is a subtle 
change compared either to the albedo change values of 0.1 used in GCM desertification experiments 
(6), or to what could be estimated from satellite observations in earlier decades. Current satellite 
systems will be capable of measuring this level of variation for future decadal fluctuations. The 
present results suggest the importance of such measurements, due to the role vegetation feedbacks 
can play in interannual/interdecadal climate variability in climatically sensitive zones like the Sahel. 
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List of Figures 

1 Annual rainfall anomaly in mm y -1 (vertical bars) over the West African Sahel ( 13N- 
20N, 15W-20E) from from 1950 to 1998: (A) observations from Hulme (1) (B) model 
with non-interactive land-surface hydrology (fixed soil moisture) and non-interactive 
vegetation (SST influence only, LoffVoff); (C) model with interactive soil moisture 
but non-interactive vegetation (LonVofF); (D) model with interactive soil moisture 
and vegetation (Lon Von); showing the role of interactive soil moisture and vegetation 
in enhancing the SST-induced interdecadal rainfall variability from the 1950s to the 
1990s, and also modifying the interannual variability. Also plotted (as connected 
circles; labeled on the right) are: (A) the Normalized Difference Vegetation Index 
(NDVI) (28); (C) model simulated annual soil moisture anomaly (mm); (D) model 
simulated leaf area index (LAI) anomaly. All the anomalies are computed relative 

to the 1950-98 base period 11 

2 Sahel rainfall difference between the period 1950-67 and the period 1968-87 for 
observations (filled square) and three 9-memeber ensemble runs with and without 
interactive soil moisture and vegetation, similar to the ones in Fig. 1B-D. Open 
circles denote individual ensemble members with different initial conditions and 
hence different chaotic internal variability. Crosshairs denote ensemble means. In 

mm y” 1 12 

3 Response of rainfall to a sinusoidal forcing in the idealized linear system (Eqs. 4, 

5), illustrating the dependence of vegetatioin feedback on forcing frequency. The 
response amplitude (P / / Fo, solid line) and phase lag (^>, in radians, dashed fine) are 
plotted as a function of the forcing frequency ojt (normalized by the vegetation time 
scale r) 
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Figure 1: Annual rainfall anomaly in mm y -1 (vertical bars) over the West African Sahel (13N-20N, 
15W-20E) from from 1950 to 1998: (A) observations from Hulme (1) (B) model with non-interactive 
land-surface hydrology (fixed soil moisture) and non-interactive vegetation (SST influence only, 
LofFVofF); (C) model with interactive soil moisture but non-interactive vegetation (LonVoff); (D) 
model with interactive soil moisture and vegetation (Lon Von); showing the role of interactive soil 
moisture and vegetation in enhancing the SST-induced interdecadal rainfall variability from the 
1950s to the 1990s, and also modifying the interannual variability. Also plotted (as connected 
circles; labeled on the right) are: (A) the Normalized Difference Vegetation Index (NDVI) (28); 
(C) model simulated annual soil moisture anomaly (mm); (D) model simulated leaf area index 
(LAI) anomaly. All the anomalies are computed relative to the 1950-98 base period. 
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Figure 2: Sahel rainfall difference between the period 1950-67 and the period 1968-87 for 

observations (filled square) and three 9-memeber ensemble runs with and without interactive soil 
moisture and vegetation, similar to the ones in Fig. 1B-D. Open circles denote individual ensemble 
members with different initial conditions and hence different chaotic internal variability. Crosshairs 
denote ensemble means. In mm y -1 . 
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Figure 3: Response of rainfall to a sinusoidal forcing in the idealized linear system (Eqs. 4, 5), 
illustrating the dependence of vegetatioin feedback on forcing frequency. The response amplitude 
(P'/Fq, solid line) and phase lag (<f>, in radians, dashed line) are plotted as a function of the forcing 
frequency lot (normalized by the vegetation time scale r). 



